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Abstract: Sand fences made of punched steel plate (PSP) have recently been applied to control 
wind-blown sand in desertified and Gobi areas due to their strong wind resistance and convenient i situ 
construction. However, few studies have assessed the protective effect of PSP sand fences, especially 
through field observations. This study analyzes the effects of double-row PSP sand fences on wind and 
sand resistance using field observations and a computational fluid dynamics (CFD) numerical simulation. 
The results of field observations showed that the average windproof efficiencies of the first-row and 
second-row sand fences were 79.8% and 70.8%, respectively. Moreover, the average windproof 
efficiencies of the numerical simulation behind the first-row and second-row sand fences were 89.8% and 
81.1%, respectively. The sand-resistance efficiency of the double-row PSP sand fences was 65.4%. Sand 
deposition occurred close to the first-row sand fence; however, there was relatively little sand on the 
leeward side of the second-row sand fence. The length of sand accumulation near PSP sand fences 
obtained by numerical simulation was basically consistent with that through field observations, indicating 
that field observations combined with numerical simulation can provide insight into the complex 
wind-blown sand field over PSP sand fences. This study indicates that the protection efficiency of the 
double-row PSP sand fences is sufficient for effective control of sand hazards associated with extremely 
strong wind in the Gobi areas. The output of this work is expected to improve the future application of 
PSP sand fences. 
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1 Introduction 


Researcher divided sand-control measures into solidification, resistance, dredging, and diversion 
categories according to their functions and properties (Bruno and Fransos, 2015; Li and Sherman, 
2015; Horvat et al., 2021). Among these, sand-resistant engineering measures primarily include 
high vertical sand barriers, sand fences, and forest belts (Hewes, 1981; Vercauteren et al., 2006). 
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Sand-resistant engineering is generally deployed at the front of sand protection system to increase 
resistance to wind-blown sand by reducing wind velocity and blocking, and intercepting sand on 
the windward side of sand fence (Snyder and Pinet, 1981; Sherman and Nordstrom, 1994; Miller 
et al., 2001; Grafals-Soto and Nordstrom, 2009; Jackson and Nordstrom, 2011). 

In recent years, a variety of sand-resistant materials (polyethylene (PE) mesh, nylon mesh, 
plug-in concrete walls, sleeper-style windscreens, etc.) suitable for the aeolian environment have 
been applied in sand-control engineering, especially in China (Yang et al., 2020); and much 
experience pertaining to sand control under different conditions has been attained. Specifically, 
sand fences made of PE and nylon mesh are more sensitive to wind velocity, and are generally 
used in low-wind energy areas where velocity is lower than 10 m/s (Yang et al., 2020). Wang et al. 
(1999) conducted a wind tunnel test to assess the protection efficiency of white nylon mesh and 
ultraviolet (UV)-resistant mesh sand fences on the Tarim Desert Highway, China. They revealed 
that the porosity of white nylon mesh sand fence was 64%, and wind velocity on the leeward side 
of the fence was approximately 20% lower than the initial wind velocity, while the porosity of 
UV-resistant nylon mesh sand fence was 58%. Qu et al. (2001) found that wind velocity on the 
leeward side of a sand fence was 30% lower than the initial wind velocity, and sand deposition 
efficiency of a nylon mesh sand fence with a porosity of 40%—45% was 50% under extreme wind 
velocity (24 m/s). In high wind energy areas, especially in high-altitude plateaus, sand fences 
such as plug-in concrete walls and sleeper-style windscreens are often used for railway protection 
from wind-blown sand (Bruno et al., 2018a, b). Yang et al. (2020) studied the sand control effects 
of different types of sand-resistant measures using wind tunnel tests. They found that plug-in 
concrete walls and sleeper-style windscreens are the least sensitive to wind velocity, and the 
wind-proof efficiency of sand fences is the highest under high wind velocity (Huang et al., 2019). 
However, due to the large pore spacing of plug-in concrete walls and sleeper-style windscreens, 
sand accumulation on the railway line on the lee side of the windshield wall remained serious, 
indicating the low sand-resistance efficiency of the windshield (Huang et al., 2019). In addition, 
plug-in concrete walls and sleeper-style windscreens have many disadvantages, such as 
inconvenient transport, tedious construction, and poor economic efficiency. In contrast, punched 
steel plate (PSP) sand fences proposed in recent years are made of a new type of sand-resistant 
material, and can be conveniently and quickly installed during the construction process. Wang et 
al. (2018) tested the protection provided by PSP sand fences using wind tunnel experiments, and 
found that windproof and sand-resistance efficiencies of the single-row and double-row PSP sand 
fences were greater than 70% and 80%, respectively, indicating that PSP sand fences could 
effectively control sand hazards caused by extremely strong wind in the Gobi areas. However, to 
date, only wind tunnel tests have been conducted to assess the protective effect of PSP sand 
fences, and verification through field observation data is lacking. Due to space limitations, wind 
tunnel tests hardly produce the correct geometry, kinematics, and dynamics simultaneously (Li 
and Sherman, 2015). Therefore, the results of wind and sand flow around sand fences are quite 
different from those in the field (Hagen et al., 1981). 

The Golmud—Korla Railway has a total length of 1214 km, and extends from the city of 
Golmud in Qinghai Province in the east to the city of Korla in Xinjiang Uygur Autonomous 
Region in the west. The third railway out of Xinjiang, the Golmud—Korla Railway improves the 
railway network in western China, and has great significance for promoting the development of 
the region around the line by strengthening national unity and enhancing the convenience of 
transportation from the mainland of China to central Asia and southern Europe (Zhang et al., 
2019a). The location of the railway line in the Qaidam and Tarim basins primarily on the leeward 
side of strong wind velocity with stable wind direction for long periods and seasonality. Moreover, 
the underlying surface on both sides of the railway is primarily sandy gravel, which easily 
releases wind-blown sand under strong wind. Consequently, control of wind-blown sand over the 
Gobi areas is a key technology for railway safety (Bruno et al., 2018c; Tan et al., 2020). The 
geographical and environmental characteristics along the Qinghai Province of the Golmud—Korla 


JOURNAL OF ARID LAND 


Railway make the use of plants for sand fixation unfeasible (Zhang et al., 2019b). Thus, 
sand-resistant measures are primarily utilized for railway protection, and double-row PSP sand 
fences are installed in the Zhongzaohuo area along the railway line. However, an evaluation of the 
shelter efficiency of the double-row PSP sand fences, particularly through in situ field 
observations and numerical simulation, has not been reported in the literature. Due to the 
complexity of wind field, observations can only measure the average wind velocity and sand flux 
density at limited locations, and hardly reveal the turbulence characteristics around the entire sand 
fence (Li and Sherman, 2015; Bruno et al., 2018a). Field observations combined with numerical 
simulations can provide greater insights into the complex wind field than wind tunnel test (Horvat 
et al., 2020). 

In this study, an evaluation of the protection effect of the double-row PSP sand fences laid in 
the Golmud—Korla Railway was conducted by field observation to address their windproof and 
sand-resistance efficiencies. In addition, the computational fluid dynamics (CFD) numerical 
simulation method was used to verify characteristics of the flow field on the leeside of the 
double-row PSP sand fences. Overall, this study provides favorable and theoretical support for the 
design and improvement of the railway sand-control system, which is of great significance for the 
safe operation of railway. 


2 Study area 


Figure la shows the annual sand-driving wind rose in the study area (the height of weather station 
was 6 m above the ground). Sand-driving wind velocity in the region around the Golmud—Korla 
Railway is approximately 6.0 m/s at a height of 6 m above the ground (Zhang et al., 2019a). The 
dominant wind directions in this region were WNW (west-north-west), W (west), and NW 
(north-west), accounting for 56.1% of the annual sand-driving wind frequency. The secondary 
wind directions were N (north) and NW (north-west), accounting for 12.6%. The maximum wind 
velocity was 21.2 m/s, with an instantaneous value of 23.1 m/s. The maximum wind velocity in 
2018 occurred to the NW (north-west). In 2018, the annual frequency of sand-driving wind in the 
region was 12.2%; the maximum monthly sand-driving wind frequency was 21.3% in May; and 
the minimum was 4.1% in December. The sand-driving wind frequency was greater than the 
average from February to July in 2018 (Fig. 1b). 

The annual drift potential (DP) of the region was calculated using the method proposed by 
Fryberger (1979). According to the calculation, the annual DP of the region was 452.07 VU, 
which indicates a high wind energy environment. The resultant drift potential (RDP) was 328.20 
VU; the resultant drift direction (RDD) was 105.56°; and the RDP/DP (the variation rate of DP) 
was 0.726, indicating a bimodal wind condition (Fig. 1c). 
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Fig. 1 Aeolian environment of study area. (a), annual sand-driving wind rose; (b), monthly frequency of 
sand-driving wind; (c), annual drift potential. DP, drift potential; RDP, resultant drift potential; VU, vector units; 
RDD, resultant drift direction. 


ZHANG Kai et al.: Sheltering effect of punched steel plate sand fences for controlling... 


Crescent dunes near the sand-control system were important sources of railway sand hazards in 
the region. Gobi wind-blown sand forms under the action of westerly winds, causing sand 
accumulation on the railbed. Seasonal rivers distributed across the region carry abundant 
sediment during the wet season, some of which settles to form a sand bed. Moreover, water 
erosion can destroy the gravel cover of the underlying surface. In the dry season, the wind 
velocity exceeds the sand-driving velocity, thus, railway wind-sand damage occurs easily. 

Double-row PSP sand fences are installed on the windward side of the railway subgrade in suit. 
The height of the fences is 2 m, and the porosity of sand fences is 50% provided by the railway 
designer. The distance between the two rows of fences is 40 m, and the fence length is 2 km. 


3 Methods 


3.1 Field observation 


HOBO U30 anemometers (U30-NRC, Onset, USA) were established at different heights. 
Additionally, 1-m-high portable sand traps were installed at 40 m (20.0 H, where H is the height 
of the sand fence) and 1 m (0.5 H) in the upwind direction of the first-row sand fence; at 1 m (0.5 
H), 20 m (10.0 H), and 39 m (19.5 H, 0.5 H before the second-row PSP sand fence) on the 
leeward side of the first-row sand fence; and at 1 m (0.5 H), 20 m (10.0 H), and 40 m (20.0 H) on 
the leeward side of the second-row PSP sand fence (Fig. 2). The heights of the HOBO 
anemometers were 0.2, 0.5, 1.0, 1.5, and 2.0 m; and the height of wind direction vane was 2 m. 
The acquisition frequency of wind velocity and wind direction was 0.1 Hz, and the recording 
interval was 1 min. The 1-m-high sand traps were a portable flat-mouth sand collector with an 
opening size of 2 cmx2 cm. Fifty sand collectors were used. The efficiency of the sand traps was 
between 72% and 87%. The wind-blown sand field was observed on 16 May, 2018 during the 
westerly wind transport process. The dominant wind direction in this region was westerly wind. 
The average angle between the observed wind-blown sand and the trend of system section was 
92°, and they were basically vertical. The wind velocity at a height of 2 m above the ground was 
10.3 m/s, and the instantaneous maximum wind velocity was 17.5 m/s. 
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Fig. 2 Field observation made on double-row PSP sand fences. (a), cross-section diagram of the anemometers; 
(b), cross-section diagram of the sand traps. PSP, punched steel plate, H, the height of the sand fence, which is 2 
m. 


3.2 Simulation method 


3.2.1 Mathematical models 
The simulation was conducted using the Reynolds number average Navier-Stokes equation, which 
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is suitable for simulating the time-average flow characteristics of long-term morphological 
dynamics around railways. Since the SST (shear-stress-transport) k-œ turbulence model has been 
proved to be an accurate model of the wind-sand solid barrier (Horvat et al., 2020), it is selected 
for the current study. All simulations were performed under steady-state conditions, and the inlet 
airflow was turbulent, incompressible, and separated. 

The entire governing equation expressed in Einstein symbols and Cartesian coordinates is as 
follow (Menter et al., 2003; Bruno et al., 2018c; Horvat et al., 2021): 


ðu; 
~=0, 1 
m2 (1) 
, . Ou; 
7 ou, 1 © ð ey Ou; | a (2) 
OX; pox; OX; OX; OX; 
Ok ð ôk * 
u; ax, = ax, os Hv) a. +P. -B ko, (3) 
gp Neg ei ng Or ep) me LA (4) 
10X; OX, OX ; k o OX; OX; 


where u; is the averaged velocity in the i direction (m/s); X;is the i direction along the X axis; uj is 
the averaged velocity in the j direction (m/s); X; is the j direction along the X axis; P is the 
averaged pressure (Pa); p is the air density (kg/m); v is the air kinematic viscosity (m7/s); k is the 
turbulent kinetic energy (m?/s?); œ is the specific dissipation rate (m7/s”); and v is the turbulent 
kinematic viscosity (m/s). The kinetic energy production term / is modeled by introducing a 


limiter to prevent the build-up of turbulence in stagnation regions: 

Ou; Ou ; Ou; 
+——., 

Ox; Ox; Ox; 

For the sake of conciseness, the definition of the blending function F; and the values of the model 

main constants B*, Ok, Go, Ge2, Co: and Coz are omitted herein. The indices can be found in detail 

from Menter et al. (2003). 

In order to reflect the influence of the roughness elements on the ground wind field, the SST 
k-@ model needs to supplement the sand-grain roughness wall function. This supplement is based 
on the equivalent sand grain roughness length (Ks) and roughness constant (Cs). Rough wall 
modification with roughness height is determined as K;=9.793z0/Cs, where C.=0.5 is the 
roughness constant; and zo is the aerodynamic roughness length (m). 

A 3D model was used for numerical simulation, which was built by using AutoCAD software. 
The calculation domain included the flat ground on the windward and leeward of the sand fence, 
as well as the sand fence itself. The length and height of the computational domain of the fluid 
were Ly=100.0 H and Lz= 20.0 H, respectively (Fig. 3). We set the size Ly=10.0 H according to 
the method of Horvat et al. (2021). The distances of the entire computational domain were 20.0 H 
and 40.0 H from the sand fences in the windward and leeward directions, respectively, which 
were sufficiently large (>15.0 H) to avoid influences of boundary conditions on the results. The 
geometric model for numerical simulation was established with a PSP sand fence at a height of 2 
m. The thickness of the sand fence was 2 mm, and the length was 20 m. The holes on the sand 
fences were elliptical, with the short side of 1.0 cm and the long side of 1.5 cm. Therefore, there 
were approximately 42x10° ellipses in a row of sand fence, and the ellipses were evenly 
distributed on the sand fence. The porosity of 50% was defined as the percentage of the area of 
the holes to the total area of the fence; and PSP sand fences were spaced at 40-m intervals. 

The medium type of the model was designated as "fluid", the inlet wall surface was defined as 
"inlet", the outlet wall surface was defined as "outflow", the upper wall surface and the left- and 
right-wall surfaces were defined as "symmetry", and PSP sand fences and the lower 
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Fig. 3 Schematic diagram of the model of computational domain and boundary conditions. H, the height of the 
sand fence, which is 2 m. 


wall surface of the model were defined as "wall". The inlet wind velocity profile was obtained by 
Equation 5; and the reference wind velocity at a specified height of 2 m was 10.3 m/s, which was 
equal to the measured wind data. To simplify the calculation, we judged the wind profile to be 
similar to the wind profile on the open natural ground (Wakes et al., 2010; Jackson et al., 2011). 
The aerodynamic roughness length zo and friction wind velocity were 0.87 mm and 0.497 m/s, 
respectively, which were estimated from the wind profile of the reference anemometer array using 
the method of Bauer et al. (1992). We calculated the inlet boundary conditions for k and œ 
according to the following formula for the assumption of a balanced boundary layer (Liu et al., 
2011): 


{a ene (6) 
K Zo 

k =1.5(0 vD", (7) 

@=k" (Cy IL, (8) 


where Uz is the wind velocity at height Z (m/s); u» is the friction wind velocity (shear wind 
velocity; m/s);« is the Karman constant, which is generally 0.4 (Richards and Hoxey, 1993; 
Santiago and Martin, 2005); z is the height above the ground (m); and zo is the aerodynamic 


roughness length, which generally refers to the height at which the velocity is 0 (m). The values 
0.25 


„ are omitted herein, which can be referenced from 


of the model main constants Uavg, I, 1, C 
Ansys (2019). 

3.2.2 Simulation method 

Ansys Workbench software was used to generate a grid for the governing transport equations. 
Because of the complexity of geometric models, this study used a hybrid grid scheme with a 
hexahedral grid, which was easily adapted to the geometric models (Fig. 4). In the simulation, the 
center of the adjacent wall unit was placed in the logarithmic layer and that the logarithmic layer 
was valid. The wall function requirement on dimensionless wall unit 30<n*<200 (Blocken et al., 
2007). The grid convergence index method was used for grid sensitivity analysis and discrete 
error estimation (Celik et al., 2008; Tan et al., 2016). An important integral parameter for the 
horizontal wind velocity at 2 m above the ground was studied as the simulation objective. We 
chose three meshes for analysis, and took three values of ny, namely nw =0.00125 H, 0.0015 H, 
and 0.003 H, where nw is the height of the adjacent wall cell (m). The overall linearity factor of 
coarsening and refinement was 1.56 that could ensure the normal operation of the computer and 
found the appropriate number of grids. The selected grids were 8.82x10° (grid 1), 4.91x10° (grid 
2), and 2.26x10° (grid 3). Table 1 shows the results of the discretization error of horizontal 
velocity. The numerical uncertainty of fine-grid solution of horizontal velocity was 1.21% in these 
sites. The solution error in grid 2 was within an acceptable range. 
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Fig. 4 Simulated mesh and the mesh near the hole of sand fence 


Table 1 Calculated local order accuracy (p) and fine-grid convergence indices ( Ce, ) for horizontal velocities 


Index Ay=50 m Ay=60 m Ay=100 m Ay=120 m 
Ugriat (m/s) 6.231 15.326 3.636 6.132 
Ugia (m/s) 5.672 14.676 3.521 5.643 
Ugria3 (m/s) 6.303 15.441 3.754 6.198 

p 1.850 2.090 1.110 1.790 
GCI, (%) 1.21 0.65 2.36 1.42 


Note: u is the horizontal velocity on the X-axis at a height of 2 m above the surface; Hy is the distance from the entrance. 


A steady-state solver was used for the Fluent simulation. The finite volume element method 
was used to calculate the flow field numerically, the second-order upwind discrete scheme was 
adopted for the momentum, and the first-order upwind formula was adopted for k and œ (Lima et 
al., 2020). The flow field was initialized by the value set by the inlet boundary condition. The 
convergence criterion of the scaled residual error was set to 10% for all the variables and the 
continuity equation; and the calculation was aggregated except for the continuity equation (when 
the remaining ratio of the continuity equation reached 10~*). Additional iterations were performed 
until the residual ratio of all variables, and the continuity equation did not further reduce (Yang et 
al., 2009). 

3.2.3 Wind and sand characteristics 


The incoming wind flow used in simulations can reflect real environment in the study site. The 
aerodynamic roughness length zo and friction wind velocity were 0.87 mm and 0.497 m/s, 
respectively, which were estimated from the wind profile of the reference anemometer array using 
the method reported by Bauer et al. (1992). The adopted particle size of the sand was 0.1 mm, 
which was the result measured in the laboratory after we took the sand on site (Zhang et al., 2021), 
and was in line with the particle size range of Raffaele et al. (2016). Meanwhile, the 
inhomogeneity in the horizontal direction was obvious to the airflow in the watershed. The 
solution to this problem was mainly to measure without setting up sand fences. We mainly used 
the same method as Liu et al (2011) to measure incident flow and approach flow. In this study, the 
"incident" position is the place in an empty domain where the model (at the first-row sand fence) 
is to be settled in the corresponding simulation. The "inlet" position is the place at the entrance 
(20.0 H before the first-row sand fence). The "approach" flow is measured before the disturbed 
area of the model (10.0 H before the first-row sand fence). It can be seen from Figure 1 that the 
wind speed profiles at different locations in the basin were consistent, and the k, œ, and p profiles 
were also the same (Fig. 5). 


3.3 Data processing 


First, the wind velocities at different measuring points and different heights of the sand-control 
system cross section were measured by the multi-gradient HOBO anemometers, and the 
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Fig. 5 Inlet, approach, and incident flows in empty computational domains. (a), vertical profiles of the turbulent 
kinetic energy (k); (b), vertical profiles of the specific dissipation rate (œ); (c), vertical profiles of the absolute 
pressure (p); (d), vertical profiles of the mean wind velocity (u). 


windproof efficiency of the sand-control system (i.e., the attenuation percentage of the wind 
velocity of the flow after it entered the sand-control system) was calculated. The calculation 
formula of the windproof efficiency of the sand-control system follows (Wang et al., 2018; Zhang 
et al., 2021): 
7, = oz ~az] 100%, (9) 
uoz 

where 7 is the windproof efficiency of the sand-control system, with values ranging from 0 to 1; 
uoz is the wind velocity of the peripheral control point at height Z (m/s); and wiz is the wind 
velocity of the internal measurement point at height Z (m/s). The change of wind velocity only 
refers to the change of magnitude, and has nothing to do with direction. 

Additionally, the sand flux density at different measuring points changes with the attenuation of 
the sand-control system, that is, the sand-resistance efficiency (Wang et al., 2018), which was 
calculated as follows: 


No (2, 2a) 100% (10) 


Z 

where 7g is the sand-resistance efficiency of the sand-control system; Qz is the sand flux density 
within the height of 1 m at the peripheral control point of the sand-control system (g/(m7+s)); and 
Qzı is the sand flux density within the height of 1 m at the internal measurement point of the 
sand-control system (g/(m7+s)). 
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4 Results 


4.1 Field observation 
4.1.1 Variation of windproof efficiency 


Figure 6 shows the variation of wind velocity around the double-row sand fences under field test 
conditions. It can be seen from the figure that when the airflow passed through the double-row 
PSP sand fences, the airflow velocity decreased after the first-row sand fence, and the airflow 
velocity was significantly reduced. The airflow velocity was still lower than the original airflow 
velocity before the second-row sand fence. After passing through the second-row sand fence, the 
airflow velocity decreased again. As the distance from the second-row sand fence increased, the 
wind velocity gradually recovered. 
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Fig. 6 Wind velocity variation around double-row sand fences. X/H is the ratio of the distance in the X direction 
to H; and Vx/Vuis ratio of the velocity in the X direction to the velocity of the fence height H. 


Equation 8 was used to calculate the windproof efficiency at different wind velocities. For an 
average wind velocity of 10.3 m/s, the windproof efficiencies at heights of 0.2 and 0.5 m and at a 
distance of 0.5 H behind the first-row sand fence were 92.4% and 87.9%, respectively. At 10.0 H 
from the first-row sand fence, the windproof efficiencies decreased to 81.6% and 72.8%, 
respectively, indicating recovered airflow. The windproof efficiencies at heights of 0.2 and 0.5 m 
at 0.5 H behind the second-row sand fence increased to 90.6% and 82.4%, respectively. The 
windproof efficiencies at 10.0 H and 20.0 H behind the second-row sand fence decreased to 
72.1% and 67.3%, and 49.6% and 37.4%, respectively. Therefore, a significant wind velocity 
reduction zone behind PSP sand fences reduced the carrying capacity of wind-blown sand. The 
sand particles settled to form a sand deposit, and PSP sand fences played a protective role. 

4.1.2 Variation of sand-resistance efficiency 

Figure 7a shows the relationship curves between the sand flux density at different heights and the 
sand fence, where the sand flux density is defined as the amount of sediment passing through a 
unit area in unit time (g/(m°-s)). The sand flux density decreased with increasing height before the 
first-row sand fence at 20.0 H. Nonlinear regression analysis of the experimental data showed that 
the sand flux density followed a decreasing exponential function with height. The following 
exponential decay function was used for the relationship between sand flux density and height: 


Q, =C+Ae~?”), (11) 


where A, B and C are all constants; and Z is the height from the ground (m). 

However, for wind-blown sand passing through PSP sand fences, the sand flux density no 
longer followed the declining exponential function, and the sand flux density decreased 
significantly behind the sand fence. For wind velocity on the leeward side of the sand fence that is 
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less than the sand-driving wind velocity, the transmission between the surrounding winds is very 
weak; therefore, most sand grains settle around the sand fence, thus stopping the saltation process 
of surface sand grains. 

The sand flux density was calculated to analyze the variation along the sand-control system, 
and to evaluate the sand-resistance efficiency of the sand-control system (Fig. 7b). For an average 
wind velocity of 10.3 m/s, the average sand flux density at the control point outside the 
sand-control system during observation period was 14.3 g/(m*-s). The sand flux density decreased 
to 8.8 and 9.6 g/(m?.s) at 0.5 H and 10.0 H behind the first-row sand fence, respectively. After 
passing through the second-row sand fence, the average sand flux density at 0.5 H and 10.0 H 
behind the sand fence further decreased to 6.5 and 5.6 g/(m*-s), respectively. Finally, the average 
sand flux density at 20.0 H behind the second-row sand fence was 5.0 g/(m*-s). The average sand 
flux density leeward the double-row PSP sand fences continued to decrease, especially leeward 
the second-row sand fence. 

The sand-resistance efficiency was calculated by Equation 10 to further analyze the protection 
benefit of the sand-control system (Fig. 7c). The results showed that the sand-resistance efficiency 
was negative at 0.5 H before the first-row sand fence; thus, the mass of sand grains collected at 
this location increased, and sand grains accumulated in front of the first-row sand fence. The 
sand-resistance efficiencies at 0.5 H and 10.0 H behind the first-row sand fence were 38.8% and 
33.1%, respectively. The sand-resistance efficiencies at 10.0 H and 20.0 H behind the second-row 
sand fence were 61.2% and 65.4%, respectively. Hence, the sand-control system of double-row 
sand fences had a good sand-resistance effect, and blocked most of wind-blown sand. 
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Fig. 7 (a), sand flux density around the double-row PSP sand fences for a wind velocity of 10.3 m/s; (b), 
variation of sand flux density along the sand-control system; (c), variation of sand-resistance efficiency along the 
sand-control system. X/H is the ratio of the distance in the X direction to H. 
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4.2 Numerical simulation 
4.2.1 Flow test in the field 


A deceleration zone formed in front of the first-row sand fence, and the wind velocity decreased 
as the airflow passed through the double-row PSP sand fences (Fig. 8a). The upper airflow lifted 
to form an acceleration zone as the airflow passed through the first-row sand fence. The airflow 
below the first-row sand fence was hindered by the formation of a large elliptic backflow zone 
behind the first-row sand fence, and the direction of the airflow changed. The airflow velocity 
was relatively low in the backflow zone but decreased before the second-row sand fence (Fig. 8b). 
After passing through the second-row sand fence, the upper airflow steadily increased, and the 
lower airflow formed a small backflow zone behind the second-row sand fence. The airflow 
continued advancing. After passing beyond the second-row sand fence, the wind velocity 
gradually recovered. 

The wind velocity of 10.3 m/s decreased rapidly at 5.0 H from the front of the first-row sand 
fence (Fig. 8c). After passing the first-row sand fence, a backflow zone generated behind the sand 
fence due to the irregular airflow and the wind velocity was negative (opposite the direction of 
incoming wind velocity). The length of the backflow zone reached 12.8 H, and the maximum 
velocity in the recirculation zone was 3.9 m/s. In the backflow zone, sand grains settled quickly 
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Fig. 8 (a), Overall diagram of the flow field distribution of the double-row PSP sand fences (Z/H is the ratio of 
the distance along the Z direction to H); (b), normalized diagram of horizontal velocity distribution diagram of 
double-row PSP sand fences; (c), horizontal velocity distribution diagram of the double-row PSP sand fences; (d), 
windproof efficiency of double-row PSP sand fences. X/H is the ratio of the distance in the X direction to H. 
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because of the counter-clockwise airflow and accumulated behind the first-row sand fence. The 
backflow zone was affected by the airflow passing through the holes of the first-row sand fence, 
and gradually disappeared. The velocity increased to a maximum value at approximately 2.5 H in 
front of the second-row sand fence. The maximum velocities at 0.2 and 1.5 m above the ground 
were 2.5 and 3.0 m/s, respectively, and decreased rapidly under the influence of the second-row 
sand fence. A small backflow zone formed beyond the second-row sand fence. The length of the 
backflow zone was approximately 4.9 H, and the maximum value of Vx/Vy in the recirculation 
zone was 0.06. Finally, the wind velocity gradually returned to its initial state. Within the height 
range of the double-row PSP sand fences, the wind velocity variation showed a W-shaped 
distribution (i.e., decreasing, increasing, and returning to the initial velocity along the flow 
direction). 

4.2.2 Windproof efficiency 

We used Equation 9 to calculate the windproof efficiency of sand fences (Fig. 8d). The windproof 
efficiency gradually increased on the windward side of the first-row sand fence as the airflow 
passed through the double-row PSP sand fences. The change in windproof efficiency was 
relatively complex for airflow between the first-row and second-row sand fences but remained 
above 60%. The windproof efficiency decreased significantly with increasing distance beyond 
10.0 H on the leeward side of the second-row sand fence. Sand movement is a sand transportation 
phenomenon in which more than 70% of moving sand is concentrated within 0.3 m of the ground. 
Therefore, the windproof efficiency of the double-row PSP sand fences was studied at a height of 
0.2 m above the ground. For a wind velocity of 10.3 m/s, the windproof efficiencies at 0.5 H, 10.0 
H, and 19.5 H beyond the first-row sand fence were 99.2%, 72.5%, and 97.6%, respectively; the 
efficiencies were 99.1%, 96.1%, and 48.2% at 0.5 H, 10.0 H, and 20.0 H beyond the second-row 
sand fence, respectively. 


5 Discussion 


5.1 Sand-resistance effect 


Sand fences are generally placed at the foremost edge of a protection system with the principal 
purposes of increasing the resistance to wind-blown sand, reducing wind velocity, and blocking 
and intercepting drifting sand. In the field test, the windproof and sand-resistance efficiencies at 
0.5 H before the first-row sand fence are 24.0% and 28.2%, respectively, with a difference of 
4.2%. The windproof and sand sand-resistance efficiency at 0.5 H before the second-row sand 
fence are 45.6% and 58.5%, respectively, with a difference of 12.9%. The windproof and 
sand-resistance efficiency at 10.0 H after the second-row sand fence are 61.2% and 63.9%, 
respectively, with a difference of 2.7%. The average windproof efficiencies behind the first-row 
and second-row sand fences were 79.8% and 70.8%, respectively. The maximum sand-resistance 
efficiency behind the second-row sand fence was 65.4%. 

Wang et al. (2018) conducted a wind tunnel test to assess the windproof efficiency and sand 
flux density of double-row PSP sand fences, and concluded that the maximum average windproof 
efficiency and the maximum sand-resistance efficiency behind the second-row sand fence were 
87.5% and 84.5%, respectively. The maximum windproof efficiency and the maximum 
sand-resistance efficiency results of the wind tunnel test were 16.7% and 19.1% higher than the 
field test results, respectively. 

Previous studies show that the saltation of sand grains changes significantly in time and space, 
even on a relatively flat quicksand surface (Zhang et al., 2010; Li et al., 2015). Gares et al. (1996) 
found that the transportation of sand grains changed by 25% within a 15-min interval. Jackson et 
al. (2006) found that the transportation of sand grains may exceed 1.5 times the variability within 
5 m. In a field experiment, the inherent variability in sand transport and the turbulence around the 
sand fence improved the variability of sand-resistance efficiency (Wilson and Yee, 2003). 
Changes in wind direction during sand transport further increase the uncertainty of the sand flux 
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density, and variations in the angle between the airflow and the sand fence change the shear force 
of the sand grain surface on the leeward side of the sand fence. However, studying multi-row sand 
fences in a wind tunnel is difficult, because the size of sand grains, the length of saltation, and the 
aperture size of the sand fences are limited after scaling the sand fence. The differences in the 
saltation height and length scale of sand grains in the wind tunnel versus in the field create 
additional complexity. 


5.2 Comparison between the numerical simulation and field observation 


The numerical simulation results showed that the airflow passed through PSP sand fences and that 
the air velocity substantially decreased under the influence of the fences, which formed a 
deceleration zone on the windward side and a backflow zone on the leeward side of PSP sand 
fences. Gillies et al. (2017) summarized the airflow behind a porous sand fence, which formed a 
new boundary layer at the top edge. The backflow zone developed along the wind direction, and 
the new boundary layer returned to the equilibrium condition at a distance generally much greater 
than 10.0 H from the leeward side of the sand fence. Airflow was ejected when passing through 
an acceleration zone above the leeward side of the sand fence (Wilson et al., 1990). Thus, a 
circular or upright backflow zone was formed, and the direction of the airflow close to the surface 
was opposite to that of the incoming wind velocity. 

By comparing Figures 7 and 9, it can be seen that the wind velocity gradually decreased after 
the first-row and second-row sand fences. Taking 20.0 H behind the second-row sand fence as an 
example, the wind velocities at the heights of 0.2, 0.5, 1.0, and 1.5 m from the ground were 0.50, 
0.63, 0.69, and 0.74 times of the inlet wind velocity in the field test. In the numerical simulation, 
the wind velocities at the heights of 0.2, 0.5, 1.0, and 1.5 m from the ground were 0.52, 0.49, 0.51, 
and 0.55 times of the inlet wind velocity. The average wind velocity of the field test was 0.64 
times of the inlet wind velocity in the field test, and the average wind velocity of the numerical 
simulation was 0.52 times of the inlet wind velocity in the numerical simulation, and the former 
was 1.2 times that of the latter. 

The average windproof efficiency of the numerical simulation behind the first-row and 
second-row sand fences were 89.8% and 81.1%, respectively. The average windproof efficiencies 
of the first-row and the second-row sand fences differed by 10% and 10.4% from the numerical 
simulation results, respectively. In the numerical simulation process, referring to the work of 
Horvat et al. (2021), the ground shear force (t) is proportional to the sand particle size, so the 
value of t is 0.045 Pa. t is defined as the surface shear stress by wind flow. When 7/7; is smaller 
than —1, it represents backward erosion; when it is in the range of [-l, 1], it represents 
sedimentation; and when it is larger than 1, it represents positive erosion (abbreviated as erosion). 
Figure 9 depicts the simulated deposition and erosion around the fence. It can be seen from the 
figure that the length of the sand on the windward side of the first-rwo sand fence is 1.0 H, and 
the length of the sand on the leeward side of the fence is 4.0 H. There is also sand accumulation 
between the first-row and the second-row sand fences with a length of 6.0 H, and a sand 
accumulation of 10.3 H on the leeward side of the first-row sand fence. Hotta and Horikawa 
(1991) showed that, for a porosity percentage of 50%, the leeward length of the in situ measured 
sand pile is 10.0 H (The fence heights were 0.9 and 1.2 m). The ratio of sand accumulation length 
to fence height after the first-row sand fence is consistent with the field measurement data of 


— Erosion tit >l 
— Sedimentation l<r/r<1 
— Backward erosion t/t < -1 
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Fig. 9 Deposition and erosion around the fence. X/H is the ratio of the distance in the X direction to H. q, 
surface shear stress by wind flow; tt, ground shear force. 
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Hotta and Horikawa (1991). In terms of the pattern and length of deposition, most sand grains 
pass through the first-row sand fence and pile on the leeward side. Sand deposition occurs close 
to the first-row sand fence. However, there is relatively little sand on the leeward side of the 
second-row sand fence (Fig. 10). The length of the sand on the windward side of the first-row sand 
fence in the field ranges from 1.8 H to 2.9 H, and the length of the sand on the leeward side of the 
first-row sand fence ranges from 5.7 H to 12.6 H; the length of the sand on the windward side of 
the second-row sand fence ranges from 1.2 H to 2.5 H, and the length of the sand on the leeward 
side of the first-row sand fence ranges from 3.6 H to 8.5 H. Thus, the length of the deposited sand 
in the numerical simulation is consistent with the results of the field test. Due to the influence of 
the first-row sand fence, the sand velocity drops to the threshold velocity for sand movement, 
causing most of the sand particles accumulate behind the first-row sand fence. 

This study simulated the interaction between wind-blown sand and sand fences based on the 
integration of field observation data in the Gobi areas. The velocity around the double-row PSP 
sand fences was basically consistent with the field observation data. Moreover, numerical 
simulation provided more detailed data. 


(a) 


Fig. 10 Double-row PSP sand fences on study area. (a), first-row sand fence; (b), second-row sand fence. 


6 Conclusions 


When the wind-sand flow passed through the double-row PSP sand fences, the wind velocity 
decreased and formed a vortex after the first-row and second-row sand fences, and the wind 
velocity gradually recovered to the initial wind velocity as it moved away from the sand fence. 
The double-row perforated board sand fences played a good windproof effect in the prevention 
and control of sand hazards in the Gobi areas. The windproof efficiency of the double-row sand 
fences can reach more than 80.0%. The sand transport flux no longer followed the exponential 
decreasing function, and the sand transport flux decreased significantly after the sand fences. 
Away from the sand fences, the sediment flux and height followed an exponentially decreasing 
function. The maximum sand-resistance efficiency of the double-row sand fences was 65.4%, the 
sand accumulation was mainly distributed after the first-row sand fence, and the sand 
accumulation length of the sand fence was approximately 10.0 H. In the field observation, the 
locations and number of monitoring points are limited by the characteristics of the field 
environment, number of instruments, wind velocity, and wind direction. Based on the integration 
of field system observation data, the CFD simulation method can be used to simulate the 
interaction between wind-blown sand and sand fences to obtain more detailed data in the Gobi 
areas. 
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